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We present a systematic study of interference and fc-point sampling effects in the supercell ap- 
proach to phase-coherent electron transport. We use a representative tight-binding model to show 
that interference between the repeated images is a small effect compared to the error introduced by 
using only the F-point for a supercell containing (3,3) sites in the transverse plane. An insufficient 
fc-point sampling can introduce strong but unphysical features in the transmission function which 
can be traced to the presence of van Hove singularities in the lead. We present a first-principles 
calculation of the transmission through a Pt contact which shows that the fc-point sampling is also 
important for realistic systems. 
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The growing interest for exploring the electronic con- 
duction properties of nano-structures has led to the de- 
velopment of a variety of numerical methods aimed at 
describing phase-coherent electron transport in atomic- 
scale systems from first principles. Many of these meth- 
ods are based on a supercell approach where the system of 
interest is defined inside a finite computational cell which 
is then repeated periodically in the directions perpen- 
dicular to the transport direction while the appropriate 
open boundary conditions needed to describe the cou- 
plingto the electrodes are imposed in the parallel direc- 
tion [1 BUSHEHSHIlflQ Within the supercell 
approach, one is obviously limited to consider transport 
through an infinite array of contacts, however, in prac- 
tice this limitation is usually ignored. Whether the su- 
percell approach can provide a good description of trans- 
port through a single, isolated contact depends on the de- 
gree of interference between the repeated images. As the 
transverse dimensions of the supercell are increased the 
interference is expected to decrease and the result should 
approach that of a single contact. For practical first prin- 
ciples calculations this limit is very hard to realize due to 
the computational cost associated with an enlargement 
of the supercell and therefore relatively small supercells 
containing e.g. 3x3 atoms in the transverse plan e are 
typically used as a compromise 
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A second approximation which is also commonly used, 
is to evaluate the transmission function only at the Y- 
point of the two-dimensional Brillouin zone (BZ) of the 
plane perpendicular to the transport direction. For ex- 
ample this approximation has been applied in transport 
studies of Au H and Na ^1 wires and to molecular con- 
tacts containing di-thiol benzene ^lj. Again, when the 
transverse dimensions of the supercell are increased the 
T-point approximation eventually becomes valid. While 
the role of k-point sampling for total energy calculations 
is well understood, its influence on the transmission func- 



tion has, to our knowledge, so far not been systematically 
explored. In a simple picture the transmission function is 
similar to a density of states (DOS) and as such it should 
be more sensitive to the k-point sampling than the total 
energy which is an integral of the DOS. 

In this paper we investigate the effects of k-point sam- 
pling and interference on the transmission function by 
means of two examples. The first example is based on a 
tight-binding (TB) model. The simplicity of this model 
allows for a complete separation of interference and k- 
point sampling effects and thus provides an opportunity 
to study the relative importance of these factors for the 
transmission. Our results show that even for a small 
cell containing 3x3 sites in the transverse directions, 
interference effects are not very significant compared to 
the error introduced by an insufficient k-point sampling. 
The second example is a first principles calculation for an 
atomic Pt contact demonstrating that k-point sampling 
is essential also for realistic systems. 

We consider a tight-binding model for resonant trans- 
port through an array of identical states which we refer 
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FIG. 1: Two-dimensional illustration of the TB model. A 
supercell is indicated by the black dots. The molecular states 
have on-site energy e a and are coupled symmetrically to semi- 
infinite cubic lattices via the hopping parameter t. All on-site 
energies of the lattice are zero and the hopping connecting ith 
nearest neighbors is denoted U. 
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to as "molecular states". The leads are described by 
semi-infinite cubic lattices with zero on-site energies and 
hopping parameter between zth nearest neighbors. In 
all calculations we use the parameters t\ — 1.0, t 2 = 0.4, 
<3 = 0.1 to describe the leads. The molecular states have 
on-site energy s a and are coupled to the leads via a single 
lattice site by the hopping parameter t. The number of 
lattice sites between the periodically repeated molecular 
states is N in each direction. A two-dimensional version 
of the model is illustrated in Fig. ^ 

The transmission function is calculated from the gen- 
eral formula ^3] 



T(e) = Tr[G r c (e)T L (e)G a c (e)T R (e)}, 



(1) 



where Gq" is the retarded/ advanced Green's function 
of the central region. The coupling strength is given by 
Fl/r = K^l/r - ^ L / R ), where Y, L/R is the self-energy 
due to the left /right lead. In the TB model the central 
region consists of the (infinite) array of molecular states. 

We label all TB orbitals in the model, including the 
molecular states, by integer coordinates \n, TO, I) where 
n refers to the direction of transport (ir-axis) and m, I 
refer to the perpendicular directions (y, z-axes). We say 
that an orbital \n,m,l) belongs to the central supercell 
if to, I G {1, . . . , N}. For each orbital in the central su- 
percell we form the Bloch states 
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Due to the periodicity of the 
system in the directions perpendicular to the transport, 
all matrices entering Eq. are diagonal with respect 
to k = (k y , k z ) and consequently the total transmission 
decomposes into a sum of k-dependent transmissions, 

T(k; e) = Tr [G r c (k; e)T L (k; e)G a c (k; e)T R (k; e)] . (3) 

The transmission per supercell is then evaluated as 

1 



T(e) = 



BZ Jl.BZ 



T(k;e)dk, 



(4) 



where 57bz is the area of the first Brillouin zone. In 
practice the integral is converted into a finite sum 
/dk/Oez — * Sk^k, where Wk represent the weights 
of the discrete k-points. 

In order to distinguish between interference effects due 
to the repeated supercells and electronic structure effects 
resulting from an insufficient k-point sampling we study 
the TB model in three different limits which we denote 
by A, B and C. Limit A is the ideal situation where the 
result for the transmission function has been converged 
both with respect to the transverse dimensions of the su- 
percell (JV) and the denseness of the k-point sampling. 
This is the situation where no interference between the 



repeated molecular states occur and it thus corresponds 
to transmission through a single state. In B the trans- 
verse dimensions of the supercell are fixed to N = 3 and 
the k-point sampling is converged to account properly 
for the electronic structure. Limit C is like B except that 
the l.BZ is sampled only at the T-point. 

In Fig. we show the transmission function in each 
of the three limits for the TB model with parameters 
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FIG. 2: Transmission through an array of molecular states in 
the 3d TB model with e a — and t = 2. For A the result has 
been converged both with respect to and the number of 
k-points. For B the number of k-points has been converged, 
while C represents the case of T-point only. The difference 
between A and B is thus pure interference while B and C 
differ due to an insufficient k-point sampling. 
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FIG. 3: Left: Band diagram for the cubic TB lattice along 
the line (k x ,0,0) corresponding to the transport direction. 
The unit cell contains 3x3 sites in the transverse plane. The 
degeneracy of each band is indicated. Right: DOS corre- 
sponding to the band diagram (full line) and the total DOS 
of the TB lattice (dashed line) calculated using a converged 
set of k-points. 
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e a = and t = 2. This choice for the parameter t corre- 
sponds to rather strong coupling between the molecular 
states and the leads, and is relevant for e.g. atomic metal 
contacts and chains or molecular contacts in which the 
molecule forms a covalent bond with the leads. The fully 
converged limit A is reached for N = 10 and with 25 
k-points, where all interference and k-point sampling ef- 
fects have disappeared. The similarity between curves A 
and B shows that interference effects are present but not 
very strong, even for this small supercell with N = 3. 
On the other hand a correct description of the electronic 
structure is crucial as can be seen from the large discrep- 
ancy between curves A and C. 

In order to understand the origin of the large difference 
between curves B and C, we have calculated the relevant 
lead DOS for both cases. We use a unit cell contain- 
ing 3x3 sites in the transverse plane and calculate the 
DOS of the TB lattice using respectively the T-point and 
a converged set of k-points. The result is shown in the 
right panel of Fig. [3J The T-point DOS has much more 
structure than the total DOS, and exhibits pronounced 
peaks due to van Hove singularities at the Id band edges, 
see left panel. The strong van Hove singularity is a char- 
acteristic of Id systems and does generally not occur in 
extended bulk systems. By sampling the l.BZ in the 
transverse plane, the density of van Hove singularities is 
increased and this provides an effective smearing of the 
DOS. By comparing the T-point DOS with curve C of 
Fig. we see that the dramatic dips in the transmission 
at energies —2.6 and 1.0 coincide with two van Hove sin- 
gularities. This clearly indicates that the strong variation 
in the T-point transmission is due to the artificial one di- 
mensional description of the electronic structure of the 
lead. 



As the coupling between the molecular states and the 
leads is reduced, the transmission function narrows down 
and eventually turns into a Lorentzian shaped resonance 
located close to the bare energy, e a , with a width de- 
termined by the product of i 2 and the DOS of the un- 
coupled lead at e a . This suggests that the width of the 
resonance should be sensitive to the lead DOS, and in 
particular that the T-point approximation could break 
down when the resonance lies close to a van Hove sin- 
gularity of the lead. That this is in fact the case can 
be seen in Fig. 0] which shows the transmission in the 
three limits for e a — —2.7 and t = 0.3. The curves A 
and B are almost identical meaning that interference ef- 
fects are negligible. In contrast the T-point resonance is 
shifted too far down in energy and is about four times 
too narrow as compared to the k-point sampled result. 
We conclude that k-point sampling can be of importance 
also for weakly coupled systems. We note, however, that 
the example shown in Fig. 01 is a worst case scenario due 
to the close proximity of s a with a van Hove singularity. 

As a second example we have calculated the transmis- 
sion function for a Pt contact using first principles tech- 
niques to evaluate the Hamiltonian. More details on the 
computational method can be found in Ref. The 
supercell contains 3x3 atoms in the directions perpen- 
dicular to the transport direction as sketched in the inset 
of Fig. The transmission has been calculated using 
1, 8 and 18 irreducible k-points, and from the result we 
conclude that the transmission function has converged 
to within 5% using 8 k-points. In fact the same conclu- 
sion holds for several other systems that we have studied, 
indicating that convergence of the transmission function 
can in general be achieved using 8 irreducible k-points 



A: N and k-point converged 
B: N=3, k-point converged 
C: A/=3, T-point 
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FIG. 4: Like Fig- El except that the parameters e a = —2.7 
and t = 0.3 have been used. 
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FIG. 5: (Color online) Transmission through the Pt contact 
shown in the inset. The transmission has been calculated 
using 1, 8 and 18 k-points in the transverse BZ. 
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for a system containing 3 atoms in the transverse cell. 
Returning to Fig.[S]it can be seen that the T-point trans- 
mission has too much structure and represents a rather 
poor approximation. This is very similar to what we 
found for the TB model, except that the number of un- 
physical dips/peaks in the transmission function is larger 
for the Pt contact. This is because the effective number 
of orbitals in a cross-section of the supercell, which de- 
termines the number of Id bands in the lead and thus 
the density of van Hove singularities, is larger in the first 
principles calculation than in the TB model. 

Generally, the T-point approximation is expected to be 
valid when the distance between the van Hove singulari- 
ties is comparable to the width of the features induced in 
the transmission function by the van Hove singularities. 
It is our experience, both from TB studies and first prin- 
ciples calculations, that the transmission calculated at a 
single k-point is improved if the k-point is chosen away 
from the T-point. This is due to the lower degree of de- 
generacy among the Id lead bands at a general k-point, 
which spreads out the van Hove singularities. 

In conclusion, we have investigated the importance of 
k-point sampling as well as interference among the re- 
peated images in the supercell approach to coherent elec- 
tron transport. Using a representative TB model, it was 
found that the k-point sampling is the most important 
factor of the two, and that an insufficient k-point sam- 
pling or, in particular, a T-point only sampling, can in- 
troduce spurious features in the transmission function 
which could be wrongly attributed to physical properties 
of the contact under study. The origin of these spurious 
features was traced to the presence of van Hove singular- 



ities in the lead. 
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